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ABSTRACT 

We introduce a Bayesian method for fitting faint, resolved stellar spectra in order to obtain 
simultaneous estimates of redshift and stellar-atmospheric parameters. We apply the method 
to thousands of spectra—covering 5160 — 5280 A at resolution 1Z ~ 20, 000 —that we have 
acquired with the MMT/Hectochelle fibre spectrograph for red-giant and horizontal branch 
candidates along the line of sight to the Milky Way’s dwarf spheroidal satellite in Draco. The 
observed stars subtend an area of ~ 4 deg 2 , extending ~ 3 times beyond Draco’s nominal 
‘tidal’ radius. For each spectrum we tabulate the first four moments—central value, variance, 
skewness and kurtosis—of posterior probability distribution functions representing estimates 
of the following physical parameters: line-of-sight velocity (r>i os ), effective temperature (T e g), 
surface gravity (logg) and metallicity ([Fe/H]). After rejecting low-quality measurements, 
we retain a new sample consisting of 2813 independent observations of 1565 unique stars, 
including 1879 observations for 631 stars with (as many as 13) repeat observations. Parameter 
estimates have median random errors of cr„ los =0.88 km s -1 , <77^=162 K, <7i ogff =0.37 dex 
and <T[ Fe/ / H ]=0.20 dex. Our estimates of physical parameters distinguish ~ 470 likely Draco 
members from interlopers in the Galactic foreground. 

Key words: methods: data analysis — techniques: spectroscopic — galaxies: dwarf — galax¬ 
ies: individual (Draco) — galaxies: kinematics and dynamics — Local Group 


1 INTRODUCTION 


The Milky Way’s dwarf spheroidal (dSph) satellites include the 
smallest, oldest, nearest, least luminous and least chemically- 
enriched galaxies known. These extreme properties have made 
dSphs the focus of investigations ranging from galaxy forma¬ 
tion to cosmology and particle physics. Several generations of 
imaging and spectroscopic surveys have dramatically improved 
knowledge of dSph luminosity functions, internal structure and 
kinematics, chemical abundances and star-formation histories. 
Thus far every deeper/wider survey has revealed new surprises— 
examples include the original discovery that dSph gravit ational po¬ 
tentials are dominated by dark matter (Aaronsonj 19 831). discover¬ 
ies of ‘ultrafaint’ satell ites dWillman et al.l2005l : IZucker et alj2006l : 
iBelokurov et al.|[2007l) . the disc overy of chemo-dynamically inde¬ 
pendent stellar sub-populatio ns l lTolstov et alj200 4: Battagli a et al.l 
120061: iBattaglia et al.1 l201lh . evidence for dwarf-dwarf mergers 
dAmorisco et alj 20141) , and discoveries of extremely metal-poor 
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stars dKirbv et al.ll2008l : lFrebel et alj|201dh . Many investigators are 
now exploiting the richness of existing data sets in the attempt to 
piece together a comprehensive understanding of small-scale struc¬ 
ture formation within the Galactic neighbourhood. 


The recent explosion of data in this field has been fuelled 
by rapid advances in instrumentation. For example, the availabil¬ 
ity of multi-object, medium- and high-resolution spectrographs at 
6 — 10m telescopes has increased the yield of a night’s observ¬ 
ing from a few to a few hundred spectra. As a result, the num¬ 
ber of dSph stars with precise kinematic and chemical abundance 
measurements h as g rown by more th an an order of magnitude in 
the past decade (IBa ttaglia et ~al][2 006 ; Walk er, M ateo & Olszewskil 
l2009l:lKirbv et al.ll20ld:lTollerud et all201 lLICollins et all201 ill . 


In an effort to propel this trend still further and, more specif¬ 
ically, to build samples over the wide fields necessary to study de¬ 
pendences of stellar velocity and metallicity distributions on po¬ 
sition, since 2006 we have been using the Hectochelle spectro¬ 
graph at the 6.5-m Multiple Mirror Telescope (MMT) on Mt. Hop¬ 
kins, Arizona, to conduct a large-scale stellar-spectroscopic sur¬ 
vey of northern dSphs. During this time we have also developed 
a Bayesian analysis pipeline that fits a flexible model to each spec¬ 
trum and delivers multi-dimensional, posterior probability distri- 
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bution functions (PDFs) for redshift as well as stellar-atmospheric 
parameters such as effective temperature (T e ff), surface gravity 
(log < 7 ) and metallicity ([Fe/H]) . From the PDFs we derive not 
only central values and variances corresponding to estimates of 
each parameter, but also higher-order moments (skewness and kur- 
tosis) that let us gauge quality of a given parameter estimate. 

Here our purpose is three-fold: 1) to describe the acquisition 
and reduction of our Hectochelle spectra, 2) to introduce our spec¬ 
tral model and the derivation of posterior PDFs and summary statis¬ 
tics for model parameters, and 3) to make publicly available a new 
data catalogue containing first results from our Hectochelle survey. 
As the subject of this initial data release, we choose the Draco dSph 
(central coordinates 0:2000 =17:20:12, $2000 =+57:54:55, distance 
D — 76 ± 6 kpc, absolute magnitu de My = — 8.8 ± 0.3, halflight 
radius r h = 221 ± 19 pc d Mated 1 1 9981 : lMcConnachiell2012L and 
references therein)), for which we have acquired, reduced and anal¬ 
ysed 6107 independent Hectochelle spectra for 3265 unique stars. 
From these spectra we obtain high-quality measurements for 2813 
observations of 1565 unique stars. 


2 OBSERVATIONS AND DATA REDUCTION 

All measurement s presented here come from our observatio ns 
with Hectochelle dSzentgyorgyi et anil998l : ISzentgvorgvill2006i) , a 
bench-mounted, multi-object, fibre-echelle spectrograph at the 6.5- 
m MMT on Mt. Hopkins, Arizona. Hectochelle uses 240 fibres to 
gather light from individual sources (fibre apertures are ~ 1 ") over 
a field of diameter 1 °. 


2.1 Target Selection 

The left-hand panel of Figure |T| shows a colour-magnitude di¬ 
agram (CMD) for point sources toward Draco, generated with 
g- and i-band photometr y fro m the Sloan Digital Sky Survey 
(SDSS; Data Release 9. lAhn et alj|2012l) and corrected for extinc¬ 
tion according to the dust maps of ISchlegel et al.1 J19981) . Black 
points indicate stars projected within 10 / of Draco’s centre. These 
centrally-located stars clearly trace Draco’s red-giant and horizon¬ 
tal branches (enclosed by polygons in the figure). Blue and red 
points indicate all stars—including those at larger angular separa¬ 
tions from the centre—for which we obtained Hectochelle spectra. 
In general we tried to select spectroscopic targets from red-giant 
and horizontal branches. Obvious outliers are associated with early 
(pre-2009) Hectochelle runs, for which we selected targets, pref¬ 
erentially near Draco’s centre, based on uncalibrated photometry 
taken with the Kitt Peak National Observatory’s 0.9m telescope 
(Mosaic-1 camera) and the Bok 2.3m telescope (90Prime camera). 
For later Hectochelle runs we selected targets directly from SDSS 
photometry. The right-hand panel of Figure Q] displays projected 
spatial distributions for all stars falling in the fiducial red-giant 
and/or horizontal branch boxes overlaid in the left panel (black 
points) as well as all Hectochelle-observed targets (blue and red 
points). Red points indicate stars for which the Hectochelle data in¬ 
dicate probable Draco membership, based on the criteria discussed 
in Section|7] 

2.2 Observations 

Our MMT/Hectochelle observations took place between April 
2006 — May 2011, during which we observed 33 unique fibre 


configurations for Draco. We observed nine of these fields on two 
separate occasions, and henceforth we treat these repeat observa¬ 
tions independently. We observed with Hectochelle’s ‘RV31’ order¬ 
blocking filter, which gives coverage over 5150 — 5300 A. The 
most prominent feature in this region is the Mgb triplet, with rest 
wavelengths of 5167 A, 5173 A and 5184 A; however, the region 
also includes an assortment of Fe II lines that can provide a direct 
measure of iron abundance. Since we bin the detector by factors of 
2 and 3 in spectral and spatial dimensions, respectively, our Hec¬ 
tochelle spectra have an effective dispersion of ~ 0.1 A per pixel, 
or ~ 0.3 A per resolution element, for an effective resolving power 
of 72. ~ 20, 000. In a typical fibre configuration, ~ 200 fibres 
are allocated to science targets and ~ 40 are allocated to regions 
of relatively blank sky. For each configuration. Table [272] lists cen¬ 
tral coordinates of the field, dates of observation and total exposure 
times, which are typically the sum over 3 — 5 sub-exposures. 

For a given fibre configuration we took Th-Ar arc-lamp ex¬ 
posures both before and after the sequence of science exposures, 
and we took a quartz-lamp exposure either immediately before or 
immediately after the science exposures. For the purpose of cali¬ 
bration and correction for variations in fibre throughput, we also 
took series of exposures (including Th-Ar and quartz calibrations) 
during evening and/or morning twilight. 

2.3 Data Reduction 

Hectochelle’s detector consists of two 2048 x 4096-pixel CCDs, 
each of which is read out through two amplifiers. In order to maxi¬ 
mize the fibre-packing factor on the CCDs, Hectochelle’s fibres are 
mounted in a zigzag pattern at the slit plane. As a result, adjacent 
spectra are offset in both spatial and dispersion directions on the 
detector. The entire pattern also curves along the spatial axis, de¬ 
manding care in identifying and tracing apertures before extraction. 

We use standard IRAlO routines to process the raw images, 
to extract one-dimensional spectra and to estimate the wavelength 
solution for each spectrum obtained in each science exposure. 
The entire p r ocedu re is similar to those pr eviously described by 
IWalker et al.l d2007l) and I Mateo et al.l (2008), except that here we 
also propagate the variance associated with the count in each pixel 
of the raw images. At the outset, for every science frame (i.e., the 
image obtained in an individual science exposure), we generate 
a corresponding variance frame in which the value assigned to a 
given pixel is 

Var(pix) = C(pix)G + R 2 , (1) 

where C(pix) is the count in analogue-to-digital units (ADU), 
G « 1.03 e'/ADU is the gain of the Hectochelle detector and 
R = 3 e _ is read noise. In order to propagate variances, we process 
variance frames in the same way that we process their correspond¬ 
ing science frames (see below). Where we rescale count levels in a 
given science frame—in order, e.g., to correct for variability in fi¬ 
bre throughput— we rescale variances by (the square of) the same 
factor. Where we combine spectra via addition or subtraction (e.g., 
to combine sub-exposures or to subtract sky background), we com¬ 
pute the combined variance as the sum of the variances associated 
with the pixels contributing to the sum or difference. 

1 IRAF is distributed by the National Optical Asteonomy Observatory, 
which is operated by the Association of Universities for Research in Astron¬ 
omy, Inc., under cooperative agreement with the National Science Founda¬ 
tion. 
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Figure 1 . Left: Colour-magnitude diagram for point sources toward Draco, from SDSS photometry with extinction corrections applied. Black points indicate 
stars projected within 10' of Draco’s centre. Blue and red points indicate all stars—including those at larger angular separations from the centre—for which 
we obtained Hectochelle spectra. Red points identify probable Draco members, based on physical parameters estimated from Hectochelle spectra (Section]?]. 
Right: Spatial distribution of all stars falling in red-giant and horizontal branch selection boxes shown in the left panel (black points). Blue and red points 
again identify all Hectochelle-observed targets, and red points identify probable members. The solid green ellipse marks the boundary where Draco’s surface 
brightness falls to the background level, as estimated by Irwin & Hatzidimitrioij t i 19951) . 


For a given frame we use the IRAF package CCDPROC to 
perform overscan and bias corrections independently for each am¬ 
plifier. For each of the two CCDs, we then combine arrays from the 
two amplifiers to form a continuous image. Next we use the IRAF 
package APALL to identify the locations and shapes of spectral 
apertures, and then to extract one-dimensional spectra for science, 
quartz and Th-Ar exposures and associated variance frames. We 
initially identify aperture locations and trace patterns in the rela¬ 
tively bright quartz frames. Fixing the relative locations and shapes 
of apertures according to the quartz frames, we allow the entire 
aperture pattern to shift globally in order to provide the best match 
to the spectra detected in the corresponding science frames; typical 
shifts are < 0.01 pixels; we apply exactly the same shift to define 
apertures and traces for the Th-Ar frame. For science, quartz, ThAr 
and associated variance frames, we extract one-dimensional spec¬ 
tra from each aperture by combining (adding) counts from pixels 
along the axis perpendicular to the dispersion direction. 

Next we use the extracted quartz spectra to adjust for differ¬ 
ences in fibre throughput and pixel sensitivity. We fit the quartz 
spectrum in each aperture with a 10 th -order cubic spline, then di¬ 
vide each continuum fit by the mean value over all pixels in the 
entire frame. Then we divide each science and quartz spectrum by 
the normalized fit to the quartz continuum obtained in the same 
aperture, thereby adjusting for throughput differences. Finally we 
divide the science spectra by the throughput-corrected quartzes, 
thereby correcting for differences in pixel sensitivity. 

Next we estimate wavelength solutions, A(pix). For each ex¬ 
tracted ThAr spectrum, we use the IRAF package IDENTIFY to 
fit a 9 th -order cubic spline to the centroids of between ~ 35 — 40 
identified emission lines of known wavelength. Residuals of these 
fits typically have root mean square (rms) scatter ~ 0.006 A, or 
~ 0.3 km s -1 , similar to the minimum velocity error we determine 


from observations of solar twilight spectra (Section EH >. We as¬ 
sign the same aperture-dependent wavelength solutions to the cor¬ 
responding science frames. Except for extraction from 2D to ID, 
each spectrum retains the sampling native to the detector, such that 
wavelength solutions A(pix) generally differ from one spectrum 
to another and have non-uniform AA/Apix even within the same 
spectrum. 

Next we identify and reject pixels affected by cosmic rays. We 
fit the continuum of each wavelength-calibrated science spectrum 
with a 5 th -order polynomial, iteratively rejecting pixels with values 
that either exceed the fit by more than 3.5 times the rms of resid¬ 
uals, or are smaller than the fit by more than 1.75 times the rms 
of residuals (the smaller value of the latter tolerance excludes ab¬ 
sorption features from the fit). After ten iterations, we replace the 
variances of all high outliers with large values (lO 100 ), effectively 
removing their influence on all subsequent analysis. 


2.4 Sky Subtraction 

Next we estimate and subtract sky background from all extracted, 
sensitivity-corrected and wavelength-calibrated spectra. For the ob¬ 
served spectral range (5150 — 5300 A), telluric emission is neg¬ 
ligible and the dominant source of sky background is scattered 
sunlight. Most of our spectra were obtained in near-dark condi¬ 
tions, such that sky backgrounds are typically small. Neverthe- 
l ess, f or all_spectra we follow the sky-subtraction procedure of 
IKoposov et al.ld201 li ). Specifically, for a given frame we interpolate 
each of the N B k y ~ 30 individual (one-dimensional, throughput- 
corrected) sky spectra onto a common grid with constant spacing 
AA'/Apix' = 0.01A (i.e., oversampled by a factor of ~ 10 with 
respect to the original sampling). For each discrete wavelength of 
the oversampled sky spectrum we record the median count level 
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Table 1. Log of Hectochelle Observations of Draco fields 


Field 

“2000“ 

[hh:mm:ss] 

^2000 a 

UT Date 
[dd/mm/yyyy] 

HJD 6 

[days] 

JV ex p C 

Exp. Time d 
[seconds] 

Dra-01 

17:20:24.65 

+57:53:06.9 

19/04/2006 

2453844.87 

4 

7200 

Dra-02 

17:23:50.00 

+57:52:12.0 

25/04/2006 

2453850.78 

5 

4846 

Dra-03 

17:20:24.64 

+57:53:06.9 

23/02/2007 

2454154.98 

3 

5400 

Dra-04 

17:17:37.81 

+57:46:30.2 

27/02/2007 

2454158.96 

3 

5400 




11/03/2007 

2454170.89 

4 

7200 

Dra-05 

17:20:38.90 

+57:28:04.3 

03/03/2007 

2454162.94 

3 

5400 




09/03/2007 

2454168.90 

3 

5400 

Dra-06 

17:19:23.67 

+58:28:22.4 

06/03/2007 

2454165.90 

3 

5400 

Dra-07 

17:14:24.40 

+57:28:47.6 

22/04/2007 

2454212.89 

3 

5400 

Dra-08 

17:26:35.35 

+58:15:25.6 

23/04/2007 

2454213.83 

3 

5400 

Dra-09 

17:30:06.97 

+57:38:24.2 

24/02/2008 

2454520.94 

5 

6000 

Dra-10 

17:20:25.01 

+57:53:11.7 

27/02/2008 

2454523.99 

3 

3600 

Dra-11 

17:11:57.01 

+58:18:03.1 

27/02/2008 

2454523.93 

3 

4500 

Dra-12 

17:10:09.01 

+57:29:41.0 

20/03/2009 

2454910.85 

2 

4096 




21/03/2009 

2454911.86 

3 

7200 

Dra-13 

17:20:06.72 

+57:55:32.6 

21/03/2009 

2454911.96 

2 

4800 




24/03/2009 

2454914.92 

4 

7200 

Dra-14 

17:41:60.00 

+56:00:00.0 

23/03/2009 

2454913.87 

4 

9600 

Dra-15 

17:20:15.99 

+57:55:30.0 

15/05/2010 

2455331.74 

4 

6000 

Dra-16 

17:20:15.99 

+57:55:30.0 

15/05/2010 

2455331.84 

3 

4500 

Dra-17 

17:15:23.52 

+57:55:42.0 

16/05/2010 

2455332.73 

2 

3600 




16/05/2010 

2455332.78 

3 

4500 

Dra-18 

17:20:23.58 

+58:24:30.0 

16/05/2010 

2455332.85 

3 

4500 

Dra-19 

17:20:11.57 

+57:24:54.0 

16/05/2010 

2455332.91 

3 

4500 

Dra-20 

17:20:15.99 

+57:55:30.0 

29/01/2011 

2455590.99 

3 

3600 




25/05/2011 

2455706.92 

4 

4500 

Dra-21 

17:20:15.99 

+57:55:30.0 

25/05/2011 

2455706.77 

3 

5400 

Dra-22 

17:20:15.99 

+57:55:30.0 

25/05/2011 

2455706.84 

3 

5400 

Dra-23 

17:22:55.93 

+58:02:34.0 

26/05/2011 

2455707.75 

3 

5400 

Dra-24 

17:17:39.05 

+57:53:14.0 

26/05/2011 

2455707.83 

3 

5400 

Dra-25 

17:20:46.45 

+57:17:29.3 

26/05/2011 

2455707.91 

4 

4800 

Dra-26 

17:20:00.58 

+58:37:30.6 

27/05/2011 

2455708.86 

3 

5400 

Dra-27 

17:20:15.99 

+57:55:30.0 

27/05/2011 

2455708.93 

3 

3600 




31/05/2011 

2455712.87 

3 

3600 

Dra-28 

17:20:19.72 

+57:52:21.9 

27/05/2011 

2455708.76 

4 

9600 




28/05/2011 

2455709.89 

3 

6000 

Dra-29 

17:19:02.64 

+57:22:36.0 

29/05/2011 

2455710.84 

3 

5400 

Dra-30 

17:22:06.84 

+58:27:53.3 

29/05/2011 

2455710.91 

3 

4500 

Dra-31 

17:20:17.49 

+57:54:42.0 

29/05/2011 

2455710.76 

3 

5400 

Dra-32 

17:26:54.32 

+58:29:28.6 

30/05/2011 

2455711.90 

3 

5400 




31/05/2011 

2455712.93 

3 

3600 

Dra-33 

17:14:13.51 

+57:40:52.9 

31/05/2011 

2455712.78 

4 

7200 


a central coordinates of field 

^heliocentric Julian date at beginning of first sub-exposure 

c number of sub-exposures 

^exposure time, summed over all sub-exposures 


and estimate the variance as 2.198^(2^ ^) 1 (MA D 2 ), where 
MAD is the median abs olute deviation i Kodosov et al ] 20lH: 
iRousseeuw & Crou?3ll993l) . We then smoot m the resulting spec¬ 
trum of median sky level (and associated variance) back onto the 
real, irregularly-sampled wavelength solution that is unique to a 
given target spectrum, so that we can then subtract the median sky 
spectrum from the target spectrum, pixel by pixel. 


Finally, we combine sub-exposures (which have identical 
apertures and wavelength solutions) by taking the inverse-variance- 
weighted mean at each pixel of the sky-subtracted spectra. Figure[2] 
displays examples of the resulting Hectochelle spectra for science 
targets spanning the magnitude range 18 < g < 20.5. 


3 ANALYSIS OF SPECTRA 

We fit each sky-subtracted spectrum using a flexible model that 
adapts easily to different instrumental setups and scientific goals. 


2 We use the IDL function ‘INTERPOLATE’ for the interpolation and 
smoothing steps. 


The model is built from ‘template’ spectra that can, in principle, be 
real spectra of standard stars or synthetic spectra calculated for spe¬ 
cific stellar-atmospheric models. Requirements for templates are 1) 
they must have high signal-to-noise ratios, 2 ) their resolution must 
be similar to or higher than that of the science spectra, and 3) their 
parameters of interest (e.g., stellar-atmospheric parameters) must 
be known. 


3.1 Spectral Model 

Following iKoleva et al.l d2009b and lKoposov et al] J201 lh . we con¬ 
sider a spectral model of the form 


M(\) = Pi(\)T( \ 


1 + 


+ tq os 


( 2 ) 


where c is the speed of light. The order-/ polynomial 

A _ A I A _ A « 1 ^ 

Pi(\) =p 0 +pi 


0 

1 


1 

[ As 

+ P 2 

As 


+ ...+PI 


(A-Ao) 


A s 


(3) 
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Figure 2. Sky-subtracted Hectochelle spectra (black) for probable members of Draco (left panels) and probable foreground dwarf stars (right panels) spanning magnitudes 18 < g < 20.5. Overplotted with 
coloured curves are best-fitting model spectra synthesized from the spectral library of ide Lavemv et Hil llOlA . Text lists equatorial coordinates, SDSS g-band magnitude, and our estimates of ui os , T e ff, log g. 


[Fe/H], [a/Fe], 


gives shape to the continuum and T(A[1 + Qm( - A ] +ulos ]) is a 
continuum-normalized template that is redshifted according to a 
given 1) line-of-sight velocity, vi os , and 2) order-m polynomial, 


Qm( A) = q 0 + qi 


o 

1 


o 

1 


+ g2 



+ ■•■ + qrr 


(A — Aq) 


A s 


(4) 


The latter, wavelength-dependent redshift accounts for linear and 
higher-order systematic differences between wavelength solutions 
of target and template spectra. Notice that the zero th -order term, 
go, is degenerate with ni os ; therefore we assume in all cases that 
go = 0, leaving our velocity estimates susceptible to a uniform 
error in zero point that we can estimate using spectra from stars of 
known velocity (e.g., the solar spectra taken at twilight, see Section 

rot . 


3.1.1 Templates 

Lacking Hectochelle spectra spanning a sufficient range of stan¬ 
dard stars, we generate template spectra using a library of synthetic 
spectra that has generously been provided by Young Sun Lee (pri¬ 
vate communication, 2013), and has previously been used to es¬ 
timate stellar-atmosphe ric parame ters as part of the SEGUE Stel¬ 
lar Parameter Pipeline jLee et alJl200iM ‘SSPP’ hereafter,). The 
SSPP library contains rest-frame, continuum-normalized, stellar 
spectra computed over a grid of atmospheric parameters spanning 
4000 < Tefj/K < 10000 in effective temperature (with spacing 
ATefr/K = 250), 0 < log 10 [g/(cm/s 2 )] < 5 in surface gravity 
(Alog 10 [g/(cm/s 2 )] = 0.25 dex) and —5 < [Fe/H] < +1 in 
metallicity (A[Fe/H] = 0.25 dex). While the library spectra are 
calculated for a range in the abundance ratio of a elements to iron, 
this ratio is not independent of metallicity. Rather, [a/Fe] = +0.4 
for library spectra with [Fe/H] < — 1, [a/Fe] decreases linearly as 



























































































































































































































6 Walker, Olszewski & Mateo 


metallicity increases from — 1 < [Fe/H] < 0, and then [a/Fe] = 0 
for [Fe/H] > 0. The synthetic library spectra are calculated over 
the range 3000 — 10000 A at resolution 0.01 A /pixel, which we 
degrade to 0.05 A per pixel—twice as fine as our Hectochelle 
spectra—in order to reduce computational costs. 

We denote as Lq(X, # a tm) the original library spectrum corre¬ 
sponding to the ‘vector’ of stellar-atmospheric parameters # atm = 
(T e fj, log g, [Fe/H]). In order to account for the instrumental line- 
spread function (LSF), we smooth each library spectrum using a 
Gaussian kernel, such that the smoothed version of the library spec¬ 
trum is 


Nx,t 


^ ] Tq (Ai , 0 a tm)A 1 


T(A, d a tm, /to) — 


A;-A 


N X ,temp 


(5) 


E K ' 


Ai-A 
ho 


where AT, temp * s the number of pixels in the library spectrum and 
the kernel is 


AT 


Aj — A 
ho 


exp 


1 (Ai-A) 2 ' 

2 hi _ 


(6) 


In order to reduce computational cost further during the subse¬ 
quent fitting procedure, we initially generate three different ver¬ 
sions of each library spectrum using smoothing bandwidths ho / A= 
0.05, 0.075, 0.1. These values were determined empirically to span 
the full range of bandwidths required to fit the library’s version 
of the solar spectrum (T e fF = 5750 K, log 10 [g/(cm/s 2 )] = 4.5, 
[Fe/H] = 0) to the solar twilight spectra acquired through each 
fibre during our Hectochelle runs. 

Finally, in order to let the spectral model vary continuously 
despite the library’s coarse gridding in stellar-atmospheric param¬ 
eter space and the three discrete values for the initial smoothing 
bandwidth, we apply a second smoothing kernel over the entire 
collection of smoothed library spectra. Specifically, for any choice 
of stellar-atmospheric parameters, 0 a t m , and smoothing bandwidth, 
ho, we obtain a unique template. 


n a) = 


\ ^ t / \ n 7 \ ts- (^atmj ^atm ^0,- ^0 \ 

> v L(X, flatm* , ho. )K 2 [ - b -• - 

,-—i V h&tm h / 


Nl 

X>A 

i= 1 


^ a tm, ^atm /+, ho 


(7) 

where Nl = 25275 is the number of library spectra (including 
three smoothed versions for each of the 8425 original library spec¬ 
tra) and the kernel is 


AT 


-6> a 


ho, — hr 


itatm 

exp 


1 ({T eBi - T eff f + (log pi - logp ) 2 


/i?, 

1 «ff 


h 2 

n iogg 


([Fe/H], - [Fe/H ]) 2 + (h 0i - h 0 ) 2 ' 1 


h 2 

ft [Fe/H] 


h 2 

n h 


( 8 ) 


We set smoothing bandwidths equal to the grid spacing in each 
dimension: Ar eff = 250 K, h\ ogg = 0.25 dex, /ipe/H] = 0.25 
dex and hh = 0.025 A. Compared with alternative interpolation 
schemes that we tried, this smoothing procedure gives posterior 
probability distributions that are more Gaussian and less likely to 
cluster near the library’s grid points. 


3.2 Likelihood Function and Free Parameters 

Given the spectral model, M (A), we assume that the observed, sky- 
subtracted spectrum, S/A), has likelihood 

£(S(A)|0,si,s 2 ) = 

^_1_^[_ 1 (S(Ai) - M{\j)) 2 ' 

\/27r(siVar[S(Ai)] + s%) P . 2 siVar[S(Ai)] + ’ 

(9 

where we have introduced two additional nuisance parameters, si 
and s 2 , that rescale and add in quadrature, respectively, to observa¬ 
tional errors in order to account for cases of mis-estimated noise. 

In practice, the value of M (Ai) that we use in Equation[9]is the lin¬ 
ear interpolation, at observed wavelength Ai, of the discrete model 
we calculate from Equation [2] This interpolation is necessary be¬ 
cause a given template spectrum, T{ A), retains the discrete wave¬ 
length sampling of the synthetic library, which generally differs 
from those of the observed spectra. 

We fit to the Hectochelle spectra over the region 5160 < 
(A/A) < 5280, which gives N\ ~ 1150. For the templates we 
consider rest-frame wavelengths 5150 < A/A < 5290, giving 
AT, temp = 1400 and allowing redshifts corresponding to the range 
-575 < vios/( km s’ 1 ) < +575. 

In order to define the polynomials in Equations [3] and [4] we 
choose 1 = 5 and m = 2, respectively. These choices give suffi¬ 
cient flexibility to fit the continuum shape and to apply low-order 
corrections to the wavelength solution. We adopt scale parameters 
Ao = 5220 A and A s = 60 A, such that —1 < (A — Ao)/A s < +1 
over the entire range considered in the fits. 

With these choices the spectral model, M (A), is fully specified 
by a vector of 13 free parameters: 

0 = (thos, Teff, logs, [Fe/H], 

PO,Pl,P2,P3,P4,P5,qi,q2,ho) ■ (10) 

The first four have physical meaning and the rest are nuisance pa¬ 
rameters that specify coefficients of polynomials, smoothing band¬ 
width and adjustments to observational errors. Table 13.21 lists all 
parameters and identifies the adopted priors, all of which are uni¬ 
form over a specified range of values and zero outside that range. 


3.3 Parameter Estimation 


From Bayes’ theorem, given the observed spectrum, //(A), the 
model has posterior PDF 


p(9,si,s 2 \S{\)) 


C(S(X}\6, si, s 2 )p(9, si,s 2 ) 
p(S{ A)) 


( 11 ) 


where C(S(X)\9, si, s 2 ) is the likelihood from Equation [9] 
p(9, si, s 2 ) is the prior and 


p(S( A)) = J C(S(X)\9, si, s 2 )p{9, si, S2)d9dsrds2 ( 12 ) 

is the marginal likelihood, or ‘evidence’. 

In order to evaluate the posterior PDF we must scan the 15- 
dimensional paramete r space. For this task we use the softwar e 
package MultiNesj^ dFeroz & Hobsonl 120081 ; iFeroz et al] [2009). 


3 available at ccpforge.cse.rl.ac.uk/gf/project/multinest 
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Table 2. Free parameters and priors 


parameter 

prior 

description 

«l03/( km / S ) 

uniform between —500, +500 

line-of-sight velocity 

Teff /K 

uniform between 4000, 8000 

effective temperature 

lo sio[s7( cm / s2 )] 

uniform between 0, 5 

surface gravity 

[Fe/H] 

uniform between — 5, +1 

iron abundance 

PO 

uniform between 0, — max [S'(A)], + max[S(A)] 

polynomial coefficient (continuum; eqfl] 

Pi 

uniform between — max[S(A)], + max[S(A)] 

polynomial coefficient (continuum; eq[3l 

Pi 

uniform between — max[S(A)], + max[S(A)] 

polynomial coefficient (continuum; eq[3l 

P3 

uniform between — max[S(A)], + max[S(A)] 

polynomial coefficient (continuum; eq[3] 

PA 

uniform between — max[S(A)], + max[S(A)] 

polynomial coefficient (continuum; eq|3] 

P5 

uniform between — max[S(A)], + max[S(A)] 

polynomial coefficient (continuum; eqfl] 

ijl/(km/s) 

uniform between —10, +10 

polynomial coefficient (wavelength solution; eq.[4l 

Q2 / (km/s) 

uniform between —10, +10 

polynomial coefficient (wavelength solution; eq.Pfl 

h 0 /k] 

uniform between 0.05, 0.10 

broadens line-spread function (eq.|5] 

l°glO S 1 

uniform between —1, +6 

rescales observational errors (eq.M 

loglO s 2 

uniform between —2, +2 

adds to observational errors (eq.|9l 


a max[S(A)] is the maximum value (discounting pixels flagged as cosmic rays) of the sky-subtracted spectrum. 


MultiNest implements a nested-sampling Monte Carlo algorithm 
dSkill indl2004l) that is designed to calculate the evidence (Equa¬ 
tion EU and simultaneously to sample the posterior PDF (Equa¬ 
tion ED. Briefly, the algorithm proceeds by evolving a set of (in 
our case, 2000 ) ‘live points’ that initially are distributed through¬ 
out the parameter volume according to the prior. Then, in each 
iteration, the point with lowest likelihood is replaced by a new 
point that is drawn randomly from the prior, subject to the con¬ 
dition that the new point must have likelihood larger than that of 
the point being replaced. As the number of iterations increases, 
the set of live points becomes confined within regions of higher 
likelihood and smaller volume, allowing the evidence to be com¬ 
puted with increasing accuracy. The algorithm terminates when 
the accuracy reaches a specified threshold (we adopt a tolerance 
of 0.5 in log-evidence, with sampling efficien cy of 0 . 8, as recom¬ 
me nded by Mu ltiNest’s README file). iFeroz & Hobsonl (20081) 
and iFerozet al. (2009) demonstrate that MultiNest performs well 
even when the posterior is multi-modal and has strong curving de¬ 
generacies, circumstances that can present problems for standard 
Markov-Chain Monte Carlo (MCMC) techniques. In our own ex¬ 
periments, MultiNest consistently required fewer likelihood eval- 
uations than did the widely-used Met ropolis-Hastings algorithm 
(Metropolis et al.ll 19531 : IHastingsil 19700 . which was more prone to 
becoming trapped in the local minima generated by the presence 
of multiple absorption features. In order to fit a single Hectochelle 
spectrum, we typically require a few x 10 s evaluations of Equation 

m 


4 RESULTS 

Plotted over each example Hectochelle spectrum in Figure[2]is the 
best-fitting model spectrum, as calculated from the most likely (ac¬ 
cording to Equation |D point from the posterior sample returned by 
MultiNest. Text within each panel indicates estimates of physical 
parameters ui os (transformed to the solar rest frame), T e g, log g and 
[Fe/H]. Based on these estimates, spectra in left-hand panels corre¬ 
spond to probable members of Draco, with velocities near Draco’s 
mean of ~ —290 km s _1 , low surface gravities (log g < 3.5) char¬ 
acteristic of red giants, and low metallicities ([Fe/H] < — 2). Spec¬ 
tra in right-hand panels correspond to probable contaminants in the 
Galactic foreground, which are likely to be G dwarfs with stronger 
surface gravities (log<?>4) and metallicities ([Fe/H] > — 1). Vi¬ 
sually the foreground dwarfs have stronger absorption lines than 


do the bona fide red giants, and our spectral modelling is clearly 
sensitive to these differences. 

For the four spectra corresponding to probable Draco mem¬ 
bers in Figure [2] Figure [3] shows samples from the posterior PDF, 
as returned by MultiNest. These plots reveal various degeneracies 
among model parameters. Of the five physical parameters, T e fj, 
log <7 and [Fe/H] all show correlations. These degeneracies make 
sense because a relatively weak absorption feature can result from 
relatively low metal abundance, lower opacity in the relatively dif¬ 
fuse atmosphere that exists in a weaker gravitational field and/or 
relatively high degrees of ionization in hotter atmospheres. On the 
other hand, lack of correlation between stellar-atmospheric param¬ 
eters and wios implies that redshifts estimated from full spectra are 
relatively insensitive to the spectral types of templates. Regarding 
nuisance parameters. Figure [3] reveals correlations among polyno¬ 
mial coefficients that describe the continuum. The curving degener¬ 
acy between nuisance parameters si and S 2 indicates that the model 
prefers to adjust variance spectra either by rescaling or by adding a 
constant value, but not both. 

4.1 Posterior Probability Distribution Functions 

For a given spectrum, MultiNest’s sampling of the posterior PDF 
lets us evaluate not only the degeneracies discussed above, but also 
moments of the marginalized, 1-D posterior PDF for each model 
parameter. Given N random draws from the posterior PDF, then for 
a given parameter, X, we calculate moments of the error distribu¬ 
tion as follows. The first moment is the mean, X = N _1 x 'i< 

which we take to be the central value. The second moment is the 
variance, a 2 x = (N — l )^ 1 ^2,f =1 {xi — A') 2 , which we take to 
be the square of the la credibility interval (for a Gaussian error 
distribution, the interval X ± ax encloses 68% of the integrated 
probability). The third moment is skewness, S = N _1 JA =1 \{ x i~ 
X )/ \Ja x ] 3 , which equals zero for a symmetric distribution. The 
fourth moment is kurtosis, K = N _1 [(x t — X )/ f a' x ] 4 , 
which distinguishes Gaussian distributions (K = 3) from ‘lep- 
tokurtic’ ones with sharper peaks and fatter tails (K > 3) and 
‘platykurtic’ ones with broader peaks and weaker tails (K < 3). 

For each of our 6107 independent Hectochelle spectra of 3265 
unique stars, scatterplots in Figure [4] show how moments of the 
PDFs for physical parameters depend on the median signal-to-noise 
ratio per pixel. Each point in a given panel corresponds to a single 
spectrum. As S/N grows, the PDFs tend to become more Gaussian 
(S ~ 0, K ~ 3) and have smaller variance; thus the distributions 
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Figure 3. Correlations among model parameters, from samplings of the posterior distribution function obtained in fits to the spectra of probable Draco members shown in Figurel2l Colours correspond to those 
used for best-fitting spectra in FigurefH 


shown in Figure [4] provide objective criteria for quality control. In 
the velocity panels, we identify a cluster of 2813 measurements 
that have <r„ loa < 5 km s _1 , —1 < S < +1 and 2 < K < 4. The 
relatively small variance and near-Gaussianity of these posterior 
PDFs imply tight constraints for which the 68% credibility interval 
can be specified with reasonable accuracy by Fioa ± <r„ loa • There¬ 
fore, in all subsequent analysis and in the data tables presented be¬ 
low, we include only observations for which moments of posterior 
PDFs for Dios fall within the ranges specified above (red squares in 
the velocity panels of Figure |4]>. Because some analyses of these 
data will require only velocity information, here we do not make 


further cuts based on moments of posterior PDFs for other phys¬ 
ical parameters; however, the data tables list all four moments of 
the posterior PDF for each parameter, so users can apply any such 
cuts as needed. Furthermore, all spectra and samples from posterior 
PDFs are available for download from the following web address: 
http://www.andrew.cmu.edu/user/mgwalker/hectochelle 


4.2 Accuracy, Precision and Dependence on Templates 

In order to gauge both accuracy and precision of our estimates 
of physical parameters, we examine fits—using the same proce- 
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Figure 4. Moments of marginalized, ID posterior probability distribution functions (PDFs) of model parameters versus median signal-to-noise ratio per pixel (a resolution element is ~ 3 pixels), for 6107 
independent Hectochelle observations of 3265 unique stars. Histograms show marginalized distributions for each moment as well as S/N. Red rectangles enclose regions where posterior PDFs of velocity 
estimates are near-Gaussian (—1 < skewness < +1, 2 < kurtosis< 4) and have (J vi os ^ 5 km s — ^. Observations associated with these estimates are coloured red in other panels; subsequent analysis and 
data tables refer only to these observations. 


dure described above—to 1715 high-S/N (tens of thousands of 
counts per pixel) solar spectra acquired with Hectochelle in morn¬ 
ing/evening twilight during our Draco runs. For each physical pa¬ 
rameter, scatterplots in Figure[5]show first and (square root of) sec¬ 
ond moments of the posterior PDFs. Histograms show distributions 
of the first moments, which let us quantify empirically 1 ) mean off¬ 
sets with respect to known solar values (ui OSi g = 0, T e ft,o = 5778 
K, log lo [p 0 /(cm s -2 )] = 4.44, [Fe/H ] 0 = 0) and 2) scatter. For 
stellar-atmospheric parameters the empirical scatter is consistent 
with second moments of the posterior PDFs returned by MultiNest. 
For uios the empirical scatter of ~ 0.4 km s -1 is several times 
larger than the scatter implied by the posterior PDFs, the larger 
value reflects the typical rms of residuals in our wavelength cali¬ 
brations (Section [2.3t . For each physical parameter, X, Table |4~21 
lists the mean offset, (AX) = (X — Xq), and the standard devia¬ 
tion, a-^. 

We use these values to adjust raw parameter estimates for sci¬ 
ence targets. For each parameter we subtract the empirical offset 
from the first moments obtained from the posterior PDFs. Then 
we add the empirical scatter, in quadrature, to the (square roots of) 
second moments for all Draco observations. All subsequent dis¬ 
cussion and analysis (and data tables) refer to results that have 
been adjusted in this way. For the 2813 observations that pass the 
quality-control criteria described in Section |4~T1 our estimates of 
physical parameters have median (minimum, maximum) random 
errors of cr„ lo 3 =0.88 (0.43, 4.74) km s~\ cr Teff =162 (54, 1100) 
K, o-i og£) =0.37 (0.10, 1.49) dex and cr[ F e/H]=0.20 (0.06, 0.98) 
dex. 

We note that the offsets of our estimates from solar values are 
unique to our use of the SSPP library of synthetic spectra. We have 


Table 3. Fits to twilight spectra for three synthetic libraries 



SSPP a 

C05 b 

AMBRE C 

(^los ^los,©) 

0.81 tar 1 

0.36 km s -1 

0.76 km s —1 

(T'eff — T'eff,©) 

-303 K 

-34 K 

22 K 

<l°gS - log 90 ) 

— 0.63 dex 

—0.06 dex 

0.08 dex 

([Fe/H] - [Fe/H] 0 > 

— 0.48 dex 

0.04 dex 

—0.34 dex 


0.42 tas" 1 

0.43 km s —1 

0.43 kms _1 

a T eii 

51 K 

16 K 

43 K 

CT log 9 

CT [Fe/H] 

0.09 dex 

0.05 dex 

0.03 dex 

0.00 dex 

0.10 dex 

0.03 dex 


a library of Lee et al. i 2008 ?) 
b library of Coelht^t all 120051) 
c library ol d^Hrvem^ta!Tl2012l) 


repeated our entire fitting procedure, for both science and twilight 
spectra, using the alternative libraries o f lCoelho et alJl2005t ‘C05’ 
hereafter) and lde Lavemv et all d2012l . ‘AMBRE’ hereafter). Table 
14.21 lists solar offsets obtained for each physical parameter with 
each library. Those obtained with the C05 library are negligible 
in the sense that they are smaller than the random errors indicated 
by the posterior PDFs for the twilight spectra. The AMBRE library 
gives a significant offset from the solar value only for metallicity, 
where ([Fe/H] - [Fe/H] Q ) = -0.34 dex. 

Despite seemingly-better performances of the C05 and AM¬ 
BRE libraries, we report results only from our fits with the SSPP 
library. The irregularity and coarseness (spacing of A[Fe/H] = 0.5 
dex for [Fe/H] < — 1 and A [Fe/H] = 1 dex for [Fe/H] < —3) 
of the AMBRE grid introduce relatively large and asymmetric 
smoothing errors near the parameter space most relevant to Draco. 
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The C05 library contains no models with [Fe/H] < —2.5 and the 
resulting fits to many of our Draco spectra fail to give a lower 
bound on [Fe/H]. Reassuringly, when we compare raw estimates 
obtained from higher-metallicity ([Fe/H] ~ —2.0) Draco spec¬ 
tra, for which fits with both the SSPP and C05 library provide 
lower bounds on [Fe/H], we recover the same relative offsets— 
for each physical parameter—obtained in fits to solar spectra. This 
agreement suggests that the SSPP offsets listed in Table l4.2l do not 
depend strongly on the stellar-atmospheric parameters themselves 
and thus can be applied to our Draco spectra. Indeed, after we ap¬ 
ply these offsets our measurements for Draco stars come into good 
agreement with those of previous studies (Section[6]i. 


4.3 Repeatability 

After discarding observations for which posterior PDFs for i>i os fail 
to meet the criteria described in Section l4~Tl there remain 2813 in¬ 
dependent observations of 1565 stars. These include 1879 obser¬ 
vations of 631 stars with multiple observations passing the criteria 
described in Section ldTl with as many as 13 such observations for 
a single star. In order to examine repeatability of our parameter 
estimation, we compare results from multiple, independent mea¬ 
surements of a given star. For all science targets with repeat obser¬ 
vations, top panels in Figure[6]compare estimates derived from the 
first observation to those derived in each subsequent observation 
(Figure[7]zooms in on the velocity range of Draco members). With 
few exceptions, points scatter around 1:1 relations according to the 
variances of the posterior PDFs. 

More quantitatively, bottom panels in Figure [6] show distri¬ 
butions of deviations with respect to inverse-variance-weighted 
means—( X) = J2iLi( x i/ a x i )/Y,i'=i B (l/ (7 x i )—normalized 
by the propagated error. We do not expect these distributions to 
be Gaussian, because the calculation of weighted means introduces 
correlations among deviations that correspond to observations of 
the same star. Therefore, instead of comparing the observed dis¬ 
tributions to Gaussians, we compare them to distributions (red his¬ 
tograms in Figure[6]l that we obtain from artificial data that have the 
same number of observations for each of the same number of stars, 
but with scatter introduced only by observational errors. Compar¬ 
isons between the observed and artificial distributions show excel¬ 
lent agreement for m os , T e g and log g, indicating that the variances 
of the PDFs returned by MultiNest are indeed reliable estimates of 
68 % credibility intervals. For [Fe/H], the observed distribution has 
thicker tails than the artificial one, indicating that the variances of 
the PDFs returned by MultiNest correspond to slightly underesti¬ 
mated credibility intervals. 


5 DATA TABLES 

For all 2813 individual observations that pass the quality-control 
criteria described in Section 14.11 Table [4] lists equatorial coordi¬ 
nates, heliocentric lulian date at the beginning of the first sub¬ 
exposure, median signal-to-noise ratio per pixel, and moments es¬ 
timated from posterior PDFs for physical parameters vi os (trans¬ 
formed to the heliocentric rest frame), T e g, logg, [Fe/H] and 
[a/Fe], For all model parameters, the reported interval is X±ax — 
i.e., the central value is the first moment (mean) of the posterior 
PDF, and the credible interval is given by the square root of the 
second moment (square root of the variance). Third (skewness) and 
fourth (kurtosis) moments are listed parenthetically as superscripts 
above the credible interval. Table [5] then lists one set of estimates 



Figure 7. Same as top-left panel of Figure[6j but zoomed in on the region populated by Draco 
members. 


for each of the 1565 unique stars. For stars with multiple inde¬ 
pendent measurements. Table [5] lists the weighted mean values for 
model parameters. Users that require stricter quality-control crite¬ 
ria (e.g., cuts on moments of posterior PDFs for stellar-atmospheric 
parameters) are advised to re-compute these weighted means after 
refining Table [4] 


6 COMPARISONS WITH PREVIOUS WORK 

Spectroscopic surveys includi ng more than ~ 100 Dra co stars have 
previo usly been published bv lArmandroff et all d 1995 1 ‘A95 ’ here¬ 
after). iKlevna et all 120021 ‘K02’ hereafter) and lKirbv et al.1 d 201(1 
‘K10’ hereafter). In their line-of-sight velocity survey, A95 used 
the KPNO 4-m telescope and Hydra multifibre spectrograph (spec¬ 
tral range 4720 — 5460 A, resolution ~ 1.2 A per pixel) to measure 
259 velocities, with median uncertainty ~ 4 km s -1 , for 159 red 
gi ant candidates in Draco. A95 then merged their sample with that 
of lOlszewski et all d 1 9951) . who used the ‘old0 MMT’s original 
echelle spectrograph and photon-counting reticon (spectral range 
5160 — 5213 A, resolution ~ 0.2 A per pixel) to measure 69 veloc¬ 
ities, with median uncertainty ~ 1.7 km s -1 , for 24 probable Draco 
members. K02 used the 4.3-m William Herschel Telescope and 
WYFFOS fibre spectrograph (spectral range 8200 — 8800 A, reso¬ 
lution ~ 0.6 A per pixel) to measure 159 velocities, with mean pre¬ 
cision ~ 2 km s^ 1 , for 159 probable Draco members (K02 do not 
tabulate results for an additional 43 probable non-members). K10 
used the 10-m Keck-II Telescope and DEIMOS multi-slit spectro¬ 
graph (spectral range 6400 — 9000 A, resolution ~ 1.2 A per res¬ 
olution element) to measure velocities^ and multi-element abun¬ 
dances (including [Fe/H], [Mg/Fe], [Ca/Fe] and [Ti/Fe]) for 297 
probable Draco members. 

Figure[8]compares our velocity estimates to those of A95, K02 


4 At the time, the MMT consisted of six 1.8-m mirrors, for an effective 
diameter of 4.5 m. 

5 While K10 did not publish velocity measurements, Evan Kirby (private 
communication) has made their measurements available for the purpose of 
comparing to our results. 
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Figure 5. Estimates of physical parameters for 1715 solar spectra acquired with Hectochelle in moming/evening twilight during our Draco runs. Scatterplots show first and (square root of) second moments 
of posterior probabillity distribution functions; histograms show distributions of first moments. Text indicates 1) mean offset of first moments (means) with respect to known solar values (fios,© = 0. 
Teff,© = 5778 K, log 10 [g/cm/s 2 )] = 4.44, [Fe/H] = 0), and 2) square root of second moments (variances) of the distribution of first moments. We subtract the former from first moments of PDFs 
obtained for our science spectra, and we add the latter (in quadrature) to square roots of second moments of PDFs obtained for our science spectra; these adjustments are included in the values reported in TablePI 
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Figure 6. Repeatability of parameter estimates, from 1879 independent observations of 631 unique stars. For each (physical) model parameter, scatterplots in the top row indicate estimates obtained from 
later observations against those obtained from the first observation. Straight lines indicate 1:1 relations, and errorbars in the lower-right corner of each panel indicate median (red) and maximum (black) credibility 
intervals associated with the plotted points. Histograms in the bottom row indicate distributions of deviations with respect to the weighted mean, ( X) = Y2i=i (-X’j / a x ■ )/ cr^ 2 , normalized 

by credibility intervals, ax •» for -/V 0 bs independent measurements of each star with repeat observations. Dotted red histograms indicate distributions obtained from artificial data having the same numbers of 
observations for each star as in the real data, but with ‘observed’ values scattered only according to measurement errors. 


and K10, for which our sample has in common 52, 76 and 109 stars, 
respectively. Top panels plot our measurement^ directly against 
the earlier ones; to each of these relations we fit straight lines with 
slopes of unity, effectively allowing for systematic offsets in zero 

6 For stars with multiple Hectochelle observations, we combine indepen¬ 
dent estimates using the weighted mean. 


point. With respect to the A95, K02 and K10 surveys, these fits sug¬ 
gest zero-point offsets of Aui os = vi os ,Hecto — vi os , other =0.31 ± 
0.23 km s-\ -0.01 ± 0.17 km s” 1 and 0.02 ± 0.25 km s“\ re¬ 
spectively. 

Bottom panels in Figure[8]show histograms of velocity differ¬ 
ences, (uioa.Hecto — Auios) — vi os ,other, normalized by combined 
errors, er = . / a?„ „ . + a?,, . While our median random 

7 \/ t 'los,Hecto u los,other 
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Table 4. Hectochelle data from individual observations 0 


Q! 2000 

[hh:mm:ss] 

^2000 

HJD 

[days] 

S/N b 

^los 

[kms -1 ] c 

T eff 

[K] 

lo Sl0 9 
[dex] rf 

[Fe/H] 

[dex] 

17:20:40.31 

+57:56:18.1 

2455331.74 

3.1 

-293.4 ± 0.8t -01 ’ 3 ’ 3 ) 

4806 ± 199 (0 ' 1|3 ' 3) 

0.9 ± 0.2 (1 - 2 ' 4 - 4 ) 

-1.63 ± O.25t 0 0 ’ 3 0 ) 

17:21:17.89 

+57:56:21.2 

2455331.74 

4.4 

-292.0 ± o.7( 01 . 3 0 ) 

4613 ± 152(°- 1 ' 2 ' 8 ) 

1.1 ± o.3t°- 5 ' 2 ' 9 ) 

-1.94 ± 0.18t°- 3 ’ 3 0 ) 

17:23:26.43 

+58:08:10.7 

2455331.74 

1.6 

-37.3 ± 1.2 (01>3 - 3 ) 

4779 ± 201 (0 '°’ 2 ' 7) 

4.0 ± 0.6t 01 ' 2 ' 9 ) 

1.01 ± 0.29t -0 - 4 ’ 2 ' 5 ) 



2455706.84 

7.6 

-34.5 ± 0.5t -00 ’ 3 ’ 0 ) 

5089 ± 65( 01 ’ 3 - 0 ' 

5.0 ± o.2t 01,3 - 2 ) 

0.43 ± 0.08t°- 0 ’ 31 ) 



2455707.75 

11.5 

-34.4 ± 0.5 (o '°’ 2 ' 9 > 

5124 ± 64( 0 ' 1 " 3 ' 0 ) 

4.8 ± o.lt 01 ’ 30 ) 

0.11 ± o,07t -01 ’ 3 - 0 ) 



2455710.91 

6.7 

-33.9 ± 0.5t -00 ’ 3 ’ 0 ) 

5008 ± 69< 01 ’ 31 ) 

4.8 ± 0.2t 01 ’ 3 0 ) 

0.25 ± 0.08t o o ’ 3 -°) 

17:20:27.38 

+57:56:52.4 

2455331.74 

8.6 

-294.5 ± 0.6( 00 ’ 3 ’ 0 ) 

4620 ± 106< -0 ’ 2 ’ 3 ’ 2 ) 

0.8 ± 0.2t 10 ’ 3 s ) 

-2.23 ± 0.13t 0 0 ’ 3 0 ) 

17:21:44.98 

+58:10:50.5 

2455331.74 

2.6 

-18.7 ± l.lC—°- 1 " a - s ) 

4918 ± 118< -0 - 1 ' 2 ' 9 ) 

5.0 ± 0.3t -0 ' 4 ’ 2 - 7 ) 

-0.03 ± 0.17t -01 ’ 3 ' 3 ) 



2455708.93 

6.6 

-17.1 ± 0.5 (- °- 1 ’ 2 ’ 9 ) 

4885 ± 69 (_0 ' 1,3 '°i 

5.3 ± 0.2 f— °- 2 ’ 2 - 7 ) 

-0.25 ± 0.08 ( — °' 0 ’ 3 ' 0 ) 



2455712.87 

6.8 

-17.1 ± 0.5 (o ' o ’ 31) 

4842 ± 65 (0 ' 0 ’ 2 ' 9 ^ 

4.9 ± 0.2( -0 - 0 ’ 31 ) 

-0.31 ± 0.08 ( —01 ' 3 ' 0 ) 

17:20:47.88 

+58:08:14.3 

2455706.77 

5.4 

-50.6 ± 0.5 (—0 0 ’ 3 °) 

5332 ± 83 (01 ’ 3 °) 

5.0 ± o^t 0 ' 1 ’ 31 ) 

0.08 ± 0.09t°- 0 ’ 31 ) 



2455707.75 

6.2 

-50.0 ± 0.6^°' 0,3 '°i 

5320 ± got -0 - 1 ' 2 - 9 ) 

5.2 ± o.2t -0 - 2 ’ 2 - 6 ) 

-0.25 ± O.iot -01 ’ 3 0 ) 



2455708.93 

3.0 

-50.2 ± o.7^°‘ 1,31 i 

5009 ± not 0 ' 1 ’ 3 ' 2 ) 

4.2 ± 0.3( — °- 1 ’ 31 ) 

0.04 ± 0.13 (— °- 1 ’ 3 " 1 ) 



2455712.87 

3.2 

-50.8 ± 0.7^°' 0,3 '°i 

5143 ± lost 0 - 1 ’ 3 9 ) 

4.8 ± O. 3 I 0.°, 3 .Q) 

0.19 ± 0.13t 01 ’ 31 ) 

17:21:18.40 

+58:22:53.6 

2455331.74 

0.9 

-131.6 ± i.5(-0 3 ,3.2) 

4939 ± 265 (_0 - 2 ' 2 - 5 ) 

4.7 ± 0.7t — °- 7,2 ' 9 ) 

0.90 ± 0.34 ( — °' 4 ’ 2,0 ) 



2455708.86 

9.3 

-132.5 ± 0.5^ — °' 0,3 ' 1 ) 

5101 ± 69 (01 ’ 31 ) 

5.0 ± O.2( 01 > 3 ' 3 ) 

-0.51 ± 0.08 ( — °' 1 ’ 3 ' 1 ) 

17:21:47.10 

+58:21:51.3 

2455331.74 

4.1 

-7.8 ± 0.7 (0 0,2 ' 9 ' 

4744 ± 85(-0.1,3.1) 

5.0 ± O.3 (0 0 ’ 3 0 ) 

0.17 ± O.ll< 0 0 ’ 3 °) 

17:20:18.69 

+58:08:40.7 

2455331.74 

2.8 

-6.6 ± 0.9 (- °'°’ 3 ' 2 ) 

4959 ± 128t 0 - 1,3 0 ) 

4.5 ± 0.4 (o - 1 ’ 2 ' 9 ) 

-0.03 ± 0.18 ( — °' 2 ’ 3 ' 3 ) 



2455706.84 

10.5 

-7.2 ± 0.5(“°- 1 ’ 3 ' 3 ) 

5173 ± 5g(-0 2.3.0) 

5.5 ± o.lt -0 ' 5 ’ 2 - 7 ) 

-0.25 ± 0.07t -01 ’ 31 ) 



2455708.93 

7.3 

-6.8 ± 0.5 ( °' 2,3 - 2 ) 

5088 ± 64t -0 - 1 ’ 2 - 9 ) 

5.4 ± 0.2 ( — °- 3 ’ 2 ' 5 ) 

-0.30 ± 0.08 ( —0 0 ’ 3 ' 0 ) 



2455712.87 

7.3 

-7.2 ± 0.5t -00 ’ 2,9 ) 

5150 ± 62t -0 - 3 ’ 3 0 ) 

5.5 ± o.2( -0 - 7 ’ 3 0 ) 

-0.23 ± 0.08 ( —0 0 ’ 3 ' 0 ) 


“see electronic edition for complete data table. 

6 median signal-to-noise ratio per pixel 
c line-of-sight velocity in the heliocentric rest frame 
d units of g are cm s —2 


Table 5. Hectochelle data with weighted means for stars with multiple observations 0 


“2000 

[hh:mm:ss] 

$2000 

0SDSS 

[mag] 6 

*SDSS 

[mag] b 

(■^los ) 

[kms -1 ] c 

(Teff) 

[K] 

doglO 9 ) 

[dex] d 

([Fe/H]) 

[dex] 

^obs 

17:20:40.31 

+57:56:18.1 

18.94 

17.88 

-293.4 ± 0.8 

4806 ± 199 

0.9 ± 0.2 

-1.63 ± 0.25 

1 

17:21:17.89 

+57:56:21.2 

18.86 

17.78 

-292.0 ± 0.7 

4613 ± 152 

1.1 ± 0.3 

-1.94 ± 0.18 

1 

17:23:26.43 

+58:08:10.7 

18.81 

17.76 

-34.4 ± 0.3 

5066 ± 37 

4.9 ± 0.1 

0.27 ± 0.04 

4 

17:20:27.38 

+57:56:52.4 

18.07 

16.79 

-294.5 ± 0.6 

4620 ± 106 

0.8 ± 0.2 

-2.23 ± 0.13 

1 

17:21:44.98 

+58:10:50.5 

18.65 

17.50 

-17.2 ± 0.3 

4870 ± 44 

5.0 ± 0.1 

-0.26 ± 0.05 

3 

17:20:47.88 

+58:08:14.3 

19.51 

18.62 

-50.4 ± 0.3 

5230 ± 47 

4.9 ± 0.1 

0.00 ± 0.05 

4 

17:21:18.40 

+58:22:53.6 

19.04 

18.08 

-132.4 ± 0.5 

5090 ± 67 

5.0 ± 0.2 

-0.43 ± 0.08 

2 

17:21:47.10 

+58:21:51.3 

18.10 

16.79 

-7.8 ± 0.7 

4744 ± 85 

5.0 ± 0.3 

0.17 ± 0.11 

1 

17:20:18.69 

+58:08:40.7 

18.66 

17.63 

-7.0 ± 0.3 

5126 ± 34 

5.4 ± 0.1 

-0.24 ± 0.04 

4 


“see electronic edition for complete data table. 

6 from SDSS photometry, corrected for extinction 

c line-of-sight velocity in the heliocentric rest frame 

'hinits of g are cm s —2 


velocity error is less than half those reported by any of the previous 
studies, our results stand in generally good agreement with those 
of A95 and K10, with the possible exceptions of a few outliers that 
are pl ausibly explained by intrinsic velocity variabil ity (e.g., binary 
stars, Olszewski et al. 1996; Hargreaves et al. 1996). However, ve¬ 
locity deviations with respect to the K02 catalogue tend to be larger 
than can be attributed to reported errors. We note that the K02 cat¬ 
alogue does not include repeat velocity measurements, making it 
difficult to check accuracy of their reported errors. 

Figure [9] compares our estimates of stellar-atmospheric pa¬ 
rameters to those of K10, who estimate [Fe/H] from their Keck 
spectra but determine T e g from a combination of photometry 
and spectroscopy and determine log<? from photometry alone. 
K10 have already demonstrated good agreement between their 
estimates of stellar-atmospheric parameters and those_obtained 
from high-resolution Keck/HIRES spectra IShetrone et all 1200 ll ; 
IFulbright et alll2004l ; ICohen & Huanal2009h . In turn, having ap¬ 
plied the offsets described in Section 14.21 we find good agree¬ 
ment between our estimates and those of K10. Our fits of 
straight lines with slopes of unity indicate offsets of (AT e ff) = 


(Feff.Hecto 2eff,Klo) —14 i 10 K, (A log Cj') — (logyQ [i/Hecto] 
log 10 [<?Kio] =—0.05 ± 0.02 dex (cgs units) and (A[Fe/H]) = 
([Fe/H] Hecto — [Fe/H] K10 ) = 0.06 ± 0.02 dex; each of these off¬ 
sets is smaller than the smallest random error among accepted mea¬ 
surements in our data set. Histograms in lower panels of Figure [9] 
indicate that deviations between our estimates and those of K10 are 
distributed approximately according to reported errors. 


We also find good agreement for the four stars in com¬ 
mon with published HIRES da ta (one from each o f the stud¬ 
ies by Shetroneetal 1200 ll and IFulbright et al.l 120041 . two from 
ICohen & Huandl2009 ; red points in Figure [9]>. For these stars, the 
random errors that we report for [Fe/H] are about half those of 
the HIRES studies, despite the latter enjoying higher resolution, 
higher S/N and larger spectral range. We suspect that this differ¬ 
ence arises because we fit every pixel simultaneously while the 
HIRES estimates come from independent fits to individual lines. 
In any case, the direct comparisons give no indication that either 
technique underestimates errors. Given the different analysis tech¬ 
niques and range/resolution of the K10 and HIRES spectra with re- 
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Figure 8. Comparison of Hectochelle velocities to published and/or measured values for stars observed in previous studies. Top: Hectochelle velocities against those published by | Armandroff et all 119951 left) 
and lKleyna et all |2002l middle), and against unpublished velocities generously provided by | Kirby et aljl20icj right, private communication). Solid lines indicate 1:1 relation; dashed lines indicate the (constant) 
offset about which scatter is minimized. Bottom: Histograms for velocity differences—normalized by er = yj<?y Hecto a v other’ w ^ ere 17 V,Hecto and cr y other are reported errors for the two 
measurements—after applying offsets indicated by the dashed lines in top panels (these offsets are not applied to the data reported in Tables pTsl . Dotted curves indicate Gaussian distributions with unit variance. 


spect to ours, we are encouraged by the overall level of agreement 
displayed in Figure[9] 


7 RELATIONS AMONG PHYSICAL PARAMETERS 

Scatterplots in FigureflOlshow relations among physical parameters 
for the 1565 unique stars in our sample (for stars with repeat mea¬ 
surements, the displayed points indicate inverse-variance-weighted 
means). In most panels the distribution of points is clearly bimodal, 
indicating that our spectral modelling can separate Draco mem¬ 
bers from Galactic foreground in at least three dimensions: i>i os , 
log<j, [Fe/H]. The Draco members cluster near ui os ~ —290 km 
s“\ while foreground interlopers have a broader velocity distribu¬ 
tion that peaks near i>ios ~ —50 km s -1 . As metal-poor giants, 
the Draco members also have [Fe/H] and logt? estimates that are 
systematically smaller than those of the observed foreground stars, 
which tend to be relatively metal-rich dwarfs (foreground giants are 
brighter than Draco’s red giant branch and so are not selected for 
observation). 

Red boxes in the upper-left panels of Figure QO] enclose mea¬ 
surements that are closer to the centre of the Draco population than 
to the foreground population, as determined by eye. We crudely es¬ 
timate the number of Draco members in our sample by counting 
the number of stars with measurements that lie inside all the boxes. 
We count 468 such likely members. The relatively large number 
(~ 1000) of foreground contaminants is the direct consequence of 
having observed at large Draco-centric radius, where members be¬ 
come far less numerous than foreground stars but convey important 


information about the influence of the external gravitational poten¬ 
tial generated by the Milky Way. 

The two panels in the lower-right corner of Figure [TO] com¬ 
pare Draco’s colour-magnitude diagram, from SDSS photometry 
(see also Figure Q}, to a spectroscopic analogue based on our esti¬ 
mates of T e g and log g. Since log g traces magnitude only for stars 
ascending the red giant branch, Draco's RGB is clearly separated 
from foreground in the spectroscopic version. 

Finally, the right-hand panel in Figure[T]shows the spatial dis¬ 
tribution of the crudely-identified likely members, tens of which lie 
at projected distances_beyond Draco’s_nominal ‘tidal’ radius, Rt, 
as estimated by llrwin & Hatzidimitrioul 1 19951) from fits of Kina 
(l966|) dynamical m odels to star-count data. In agreement with 
iMunoz etalJ l 2()()(r i. we find likely members out to the limits of 
our data set, up to ~ 3 Rt from Draco’s centre. 


8 SUMMARY 

We have presented a new spectroscopic data set that includes ~ 
468 Draco members distributed over an area of ~ 4 deg 2 . The ob¬ 
servations include 2813 independent measurements of line-of-sight 
velocity (median random error a v ~ 0.88 km s -1 ), effective tem¬ 
perature (crT eff ~ 162 K), surface gravity (cri ogg ~ 0.37 dex) and 
metallicity (cr[ Fe / H ] ~ 0.20 dex) for 1565 stars, including 1879 
independent measurements for 631 unique stars with repeat obser¬ 
vations. These data will be useful for investigations ranging from 
galactic dynamics to chemical evolution and binary stars. 
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Figure 9. Same as Figurels! but for comparisons to stellar-atmospheric parameters published h v i K i i bv et all i20 1 (I. black points). Red points and histograms show comparisons for three stars with previously- 
published estimates from Keck/HIRES spectra iShetrone et all200lIiFulbright et ali20()4!cohen & Huan j200fft . 


All data from this study are publicly available. In addi¬ 
tion to the information provided in Tables [4] and [5] the inter¬ 
ested reader can download all spectra and samples from poste¬ 
rior PDFs (including those that do not pass quality-control cri¬ 
teria described in Section 14. 11 from the following web address: 
http://www.andrew.cmu.edu/user/mgwalker/hectochelle 
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cal Research Consortium for the Participating Institutions of the 
SDSS-III Collaboration including the University of Arizona, the 
Brazilian Participation Group, Brookhaven National Laboratory, 
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